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ABSTRACT 

The [HD0]/[H20] ratio is a crucial parameter for probing the history of water formation. So far, it 
has been measured for only three solar type protostars and yielded different results, possibly pointing to 
a substantially different history in their formation. In the present work, we report new interferometric 
observations of the HDO 42,2 — 42,3 line for two solar type protostars, IRAS2A and IRAS4A, located in 
the NGC1333 region. In both sources, the detected HDO emission originates from a central compact 
unresolved region. Comparison with previously published interferometric observations of the H3> 8 
3i,3 — 22,o line shows that the HDO and H2O lines mostly come from the same region. A non-LTE 
LVG analysis of the HDO and H^O line emissions, combined with published observations, provides a 
[HDO]/[H 2 0] ratio of 0.3 - 8 % in IRAS2A and 0.5 - 3 % in IRAS4A. First, the water fractionation is 
lower than that of other molecules such as formaldehyde and methanol in the same sources. Second, it 
is similar to that measured in the solar type protostar prototype, IRAS16293-2422, and, surprisingly 
enough, larger than that measured in NGC1333 IRAS4B. The comparison of the measured values 
towards IRAS2A and IRAS4A with the predictions of our gas-grain model GRAINOBLE gives similar 
conclusions to those for IRAS 16293, arguing that these protostars share a similar chemical history, 
although they are located in different clouds. 

Subject headings: astrochemistry — ISM: abundances — ISM: individual objects (NGC1333-IRAS2A, 
NGC1333-IRAS4A) — ISM: molecules — stars: formation 



1. INTRODUCTION 

The formation of solar type protostars is triggered by 
the gravitational collapse of dense fragments of molec- 
ular clouds, the so-called prestellar cores. In molec- 
ular clouds and prestellar cores, the low temperature 
and interstellar UV flux promote the formation of icy 
mantles around the dust grains. These mantles are 
mai nly composed of H ? Q, CO CO ?, H 2 CO, or CH 3 OH 
(see iBoogert fc Ehrenfreund I [20041 ). The two last be- 
come abundant with the freeze-out of CO. In addition, 
the CO freeze-out coupled with the cold conditions en- 
hance the abundance of the deuterated molecules (see 
ICeccarelli et~aT1 [20071 : ICaselli et al.ll2008h . 

Several theoretical studies have shown that the molec- 
ular deuteration is very sensitive to the physical condi- 
tions. For instance, the deuteration of gaseous species 
increases with the tot al density n-g and decreases with 
the t emperature (see iMillar et al.l Il989t iRoberts et al.l 
120031 ). Similarly, the deuteration of icy species formed 
on the grain surfaces by H and D atom addition re- 
actions, like H2O, H2CO and CH3OH, depends on the 
gaseous atomic [D]/[H] ratio, which also increases with 
the density and the CO freeze-out at low tem peratures 
([Cazaux et al.ll2TJlll iTaquet et al.ll20T2ll[2013[ ). 
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In theory, therefore, the deuteration of different man- 
tle species can be used to reconstruct the history of the 
ice formation and, consequently, of the protostar (e.g. 
I Taquet et al.l l2013). In practice, unfortunately, the di- 
rect measurement of the deuteration of frozen species is 
not possible. Observations of solid HDO towards proto- 
stars onl y yielded upper lim i ts (with HDO/H2O < few 
percent , iDartois et all 120031: IParise et all I2003D . How- 
ever, one can observe these species where the icy mantles 
sublimate, for example in the hot corino regions. Since 
the timescale needed to significantly alter the deutera- 
tion after the ice sublimation is longer than the typi- 
cal age of Class p r otostars (~ 10 5 v ersus ~ 10 4 yr, 
iCharnlev etlU 119971; lAndre et all I2000D . the measured 
deuteration of the gaseous mantle species likely reflects 
that in the ices prior to the sublimation. A compari- 
son between the measured and predicted deuteration in 
inte rstellar ices is, there fore, possible. 

In lTaauet et al. (2013), we did a first study by compar- 
ing the predictions of our gas-grain GRAINOBLE model 
(|Taauet et alj|2012al ) with the observations towards the 
protostar IRAS16293-2422 (hereafter, IRAS16293). This 
source displays a very high deuterium fractionation of 
formaldehyde and met hanol (with D/H rat i os of 15 and 
40 % resp ectively, see lLoinard et al.l 120011 : IParise et akl 
2002, 200 4]) and a, lower fractionation of water (0.1-3 %, 
seelButner et alj|2007t iVastel et al .ll2010HCou tens et al.l 



120121 120131 : iPersson et al.l 120131 ). We concluded that 
the lower fractionation of water with respect to that 
of formaldehyde and methanol is likely due to a differ- 
ent epoch of formation of the three species. Water is 
predicted to be mainly produced during the molecular 
cloud phase, while most of formaldehyde and methanol is 
formed during the colder and denser prestellar phase. We 
carried out a similar study using the measured deuter- 



ation of H 2 0, H 2 CO, and CH 3 OH towards the outflow 
shock L1157-B1 and concluded that this site had a sim- 
ilar sequence for th e formation of the i ce, but in a less 
dense environment (jCodella et al.ll2012l ). 

Encouraged by these two studies, we want here to ex- 
tend the analysis to other solar type protostars with the 
goal to reconstruct the formation history of their ices 
and compare it with the two previous cases. Ultimately, 
a similar study in a large sample of solar type protostars 
will provide us with a more complete picture of how the 
environment influences the chemical composition of the 
ices and will supply strong constraints to the theory. 

Although the fractionation of formaldehyde and 
methanol h as been measured toward s sever a l sola r type 
protostars (jParise et all 120021 |200H [20051 l200l) . ob- 
servational studies of deuterated water are scarce. In 
NGC1333-IRAS4B, the non detection of the HDO line 
at 225.6 GHz yields a n upper limit to the [HD01/[H 2 Q 1 
ratio of < 6 x 10" 4 (jJorgensen fc van Dishoeckll2010T ). 
In NGC1333-IRAS2A, several HDO and H 2 lines have 
been observed with single- dish telescopes. Using the 
Herschel Space Observatory, iKristensen et al.l (2010) ob- 
served a broad outflow component for several H 2 lines, 
but could not accurately estimate the w ater abundance in 
the warm compact region. In contrast. iLiu et all (|2011l ) 
derived the HDO abundance profile in the warm and cold 
regions of the envelope. However, single-dish telescopes 
also encompass the cold envelope, and the possible out- 
flow component. Complementary interferometric obser- 
vations, with arcsecond resolutions, are needed to resolve 
the emission coming from the hot corinos, where the ices 
are sublimated and where th e deuteration likely reflects 
the ice pristine deuterati on (jJorgensen fc van Dishoeckl 
[20101: IPersson et al.l[2013h . 

In this Letter, we present interferometric IRAM 
Plateau de Bure observations of the HDO 4 2 2 -4 2 3 
line at 143 GHz towards NGC1333-IRAS2A (hereinafter 
IRAS2A) and NGC1333-IRAS4A (hereinafter IRAS4A). 
These sources are located in the Perseus complex , in the 
NGC1333 cloud, whose distance is about 220 pc (|Cernisl 
1990). They were selected because they are the two 
line brightest low-mass protostars after IRAS16293 due 
to their distance and their luminosity and because in- 
terferometric obse^vaHpnsofJI^O have been recently 
obtained by iPersson et al.l ()2012[ ) towards them. The 
[HDO]/[H 2 0] ratio derived in the present work, com- 
bined with previous observations of deuterated formalde- 
hyde and methanol, will be compared wit h the predic- 
tions of our gas-grain model GRAINOBLE (jTaquet et al.l 
2013) to reconstruct the chemical history of these two 
protostars. 

2. OBSERVATIONS AND RESULTS 

The two low-mass Class protostars IRAS2A and 
IRAS4A were observed with the IRAM Plateau de Bure 
Interferometer on 2010 August 1, August 3 and 2011 
March 10 and in the C and D configurations of the array. 
Due to the proximity to each other, the two sources were 
observed in the same track. The 4 2j2 -4 2j 3 HDO transi- 
tion at 143.727 GHz and the 2 mm continuum emission 
have been obtained simultaneously using the WIDEX 
correlator, with a 1.8 GHz bandwidth centered at 143.5 
GHz, and providing a spectral resolution of 1.95 MHz (4 
km s _1 ). Phase and amplitude were calibrated by per- 



forming regular observations of the nearby point sources 
3C454.3, 3C84, and 0333+321. The amplitude calibra- 
tion uncertainty is estimated to be ~ 20 %. 

The data calibration and imaging were performed us- 
ing the CLIC, and MAPPING packages of the GILDAS 
softwarqj. Continuum images were produced by aver- 
aging line-free channels in the WIDEX correlator before 
the Fourier transformation of the data. The coordinates 
of the sources, and the size of the synthesized beams are 
reported in Table [TJ 



TABLE 1 

Coordinates, synthesized beams, continuum fluxes, and 

sizes of the observed low-mass protostars. 



Source 


IRAS2A 


IRAS4A 


RA 


03:28:55.56 


03:29:10.45 


Dec 


31:14:37.05 


31:13:31.18 


Synthesized Beam (") 


2.16 x 1.73 (25°) 


2.18 x 1.76 (25°) 


Continuum flux (Jy) a 


0.13 


1.08 


Source size (") a 


1.75 x 1.69 


2.13 x 1.72 



Note. — "Continuum fluxes and sizes are obtained from ellip- 
tical Gaussian fits in the (u,v) plane (i.e., deconvolved full width 
at half- maximum (FWHM) size). 

Figure [1] shows the maps of the continuum emission 
at 2 mm of IRAS2A and IRAS4A obtained after natural 
weighted cleaning. Parameters of the continuum emis- 
sion (flux and deconvolved FWHM size), obtained from 
elliptical Gaussian fits, are given in Table [1] For the two 
sources, since the FWHM size of the continuum emission 
is very similar to the size of the synthesized beam, the 
continuum emission is, therefore, not resolved. In par- 
ticular, IRAS4A is known to be a bin ary system with 
a 1.8" separation (jLoonev et al.l [2000h . as depicted in 
Figure [1] Although the continuum emission of IRAS4A 
is peaked at the southeast (SE) position rather than at 
the northwest (NW) position, we cannot resolve the two 
sources. 

Figure [1] shows the maps of the integrated HDO 4 22 - 
4 2j 3 line towards the two sources obtained after natural 
weighted cleaning. For both sources, the FWHM size of 
the HDO line is very similar to the size of the synthe- 
sized beam, as shown in Figure [TJ The emission clearly 
originates in compact regions, that are confined within 
the synthesized beam of the telescope. In particular, al- 
though the SE position of IRAS4A is the brightest in the 
continuum, the HDO line emission comes from the NW 
position. The spectra of the HDO transition integrated 
within the FWHM size are shown in Figure [TJ assuming 
Vlsr — 7 km s _1 for the two sources. Table [3] gives 
the flux and the brightness temperature of the HDO line 
transition integrated inside the FWHM size. 

Compari son with maps of th e H^O 3i3-2 20 line tran- 
sition by IPersson et al.l ([2012D shows that, although our 
beam is two times larger, most of the H DO and H 2 Q 
emissi ons originate in the same region. IPersson et al.l 
(|2012f) estimated the FWHM size of the H^ 8 emission, 
from a Gaussian fit in the (u, v) plane, and they found 
that most of the emission originates in a 0.8" ellipse, 
similar to t heir synthesiz e d bea m. Observational data 
derived by IPersson et al.l (120 12-) are given in Table [2j 
The HDO lines observed in this work are broader, by a 

6 The GILDAS package is available at 

|http://www.iram.fr/I RAMFR/GILDASl 
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Fig. 1. — Maps and spectra towards IRAS2A (left) and IRAS4A (right). Upper panels: HDO spectra integrated over the emission region 
of each source (middle panels). The velocity resolution is 4 km s -1 and the dashed blue lines mark the Vlsr a t 7 km s — 1 . Middle panels: 
Continuum maps at 143 GHz of IRAS2A (rms of 1.7 mjy/beam, contour levels are in steps of 4 a), and of IRAS4A (rms of 9.6 mjy/beam, 
contour levels are in steps of 4 a). Green contours depict the deconvolved full width at half-maximum size. Bottom panels: Maps of the 
HDO 42,2-42,3 line towards IRAS2A (with a rms of 2.9 mjy/beam km s _1 , contour levels are in steps of 3 <x) and in IRAS4A (with a rms 
of 2.8 mjy/beam km s — 1 , contour levels are in steps of 3 <r). Green contour s depic t the region showing a flux at half maximum of the line 
peak. The red crosses mark the source positions measured by Looncv ct al. (2000). The bottom-left ellipses represent the beam sizes. 



TABLE 2 
Parameters of the HDO and HJ, 8 lines observed towards IRAS2A and IRAS4A. 



Transition 


Frequency 


t^wp 


Mj 


Flux 




(J'T B <fo) oos 


Beam 


Telescope 


Av 


Ref. 




(GHz) 


(K) 


(s- 1 ) 


(Jy km s~ 


- 1 ) 


(K km s" 1 ) 


(") 




(km s _1 ) 






143.727 


319.2 








NG01333 IRAS2A 






HDO 4 2 ,2-4 2 ,3 


3.5 x 10" D 


0.43 




6.8 ±1.4 


2.2 x 1.8 


IRAM PdBi 


7 


1 


HDO li,o-li,i 


80.578 


46.8 


1.3 x 10~ 6 


0.37 




0.07 ±0.02 


31.2 


IRAM 30m 


3.9 


2 


HDO 2i.i-2i, 2 


241.561 


95.3 


1.2 x lO" 5 


2.2 




0.43 ± 0.05 


10.4 


IRAM 30m 


1.1 


2 


HDO 3i, 2 -2 2 ,i 


225.896 


167.7 


1.3 x 10~ 5 


2.6 




0.50 ±0.03 


11.1 


IRAM 30m 


1.2 


2 


Hi 8 O3i,3-2 2 , 


203.408 


203.7 


4.8 x 10~ 6 


0.98 




46 ±9 


0.9x0.7 


IRAM PdBi 


4.0 


3 




143.727 


319.2 








NGC1333 1RAS4A 






HDO 4 2 ,2-4 2 ,3 


3.5 x 10-° 


0.21 




3.2 ±0.6 


2.2 x 1.8 


IRAM PdBi 


6 


1 


Hi 8 3i, 3 -2 2 ,o 


203.408 


203.7 


4.8 x 10~ 6 


0.27 




13 ±3 


0.9x0.7 


IRAM PdBi 


2.9 


3 




This work, 2: 










2012). 
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Note. — The flux uncertainties include the calibration uncertainties, 


estimated to be ~ 20 %. 









factor of 2, than the other HDO and H5, 8 lines, due 
to the low spectral resolution of our observations. In 
IRAS2A, even if the bulk of the H5> 8 emission is associ- 
ated with the central warm envelope, an outflow compo- 
ne nt towards the sout hwest is also observed (see Fig. 2 
of IPersson et al.ll2012|) whereas the HDO line only orig- 
inates in the compact region located within the synthe- 
sized beam. Our maps are, therefore, in good agree- 
ment with previous single-d ish observations of H^O and 
HDO towards IR AS2A by iKristensen et all (|2010( ) and 
iLiu et all (|2011| ). described in the Introduction, which 
show that H2O mostly traces the outflow whereas HDO 
only traces the central envelope. 

In IRAS4A, both the HDO and H^O emission come 
from the NW position and the FWHM sizes are simi- 
lar to the respective synthesized beams, suggesting that 
the emission originates in the central region. Therefore, 
it is meaningful to compare the flux of the HDO and 
H^O lines originating in the FWHM central regions. 
The fluxes are given in Table [2] and are used in the next 
section to estimate the [HDO]/[H2 0] ratio. 

3. DEUTERIUM FRACTIONATION OF WATER 

3.1. Method 

A single transition line of HDO and H3, 8 does not al- 
low us to derive an accurate estimate of the HDO and 
H2O column densities towards the two protostars and, 
therefore, of [HDO]/[H2 0]. In order to derive the phys- 
ical conditions of the line emitting gas and the relevant 
column densities, w e compared the predi ctions from a 
non-LTE LVG code (|Ceccarelli et al II2003D wit h our ob- 
servations and the observations by ILiu et all f|2011h of 
several HDO lines towards IRAS2A, obtained with the 
IRAM 30m, JCMT, and APEX tel escopes. We c onsid - 
ered the collisional coefficients from iDaniel et all (J2011I ) 
for H^O and from lFaure etTaTl (|2f5l2T) for HDO. The Ein- 
stein coefficients are from the Jet Prop ulsion Laboratory 
molecular database fPic kett et al.l [19981) . 

We ran a grid of models covering a large parameter 
space in kinetic temperature T kin (15 values from 70 to 
220 K), n H (15 values from 1 x 10 6 to 1 x 10 9 cm" 3 ), 
HDO column density V(HDO) (15 values from 8 x 10 14 
to 1 x 10 17 cm -2 ), and source size 9 S (30 values from 0.1 
to 200 arcsec). In addition, we considered three values 
for the ortho to para ratio (opr) of H2: 10~ 2 (namely all 
H2 molecules are in the para state), 1, and 3 (thermal 
equilibrium value at T > 50 K). To find the best fit 
to the da t a, we excluded the 464 GHz line observed by 
ILiu et ail 1)20111 ) as it may be contaminated by the cold 
envelope emission, given its low energy level (22 K). 

3.2. Results 

IRAS2A — We ran the grid of models to reproduce the 
emission of the HDO lines towards IRAS2A. The H2 opr 
has a low influence on the column densities derived from 
the observations. Varying the H2 opr between 0.1 and 
3 causes a small variation of the results, by no more 
than 20 %, namely within the uncertainties of the ob- 
servations. In the following, we consider an H2 opr of 
3. The flux of all the HDO lines are well reproduced 
(reduced X 2 < 1) for T km ~ 75-80 K, 9 S = 0.4", and 
a wide range of «h between 6 x 10 5 and 2 x 10 8 cm~ 3 . 
The derived TV (HDO), varies between 5 x 10 17 and 10 19 



cm~ 2 and decreases with %. To evaluate the H2O col- 
umn density iV(H20) from the 203.4 GHz transition, we 
considered three physical cases that reproduce the emis- 
sion of the HDO lines (see Table [3]) . The density used 
in Case 1 (6 x 1 5 cm " 3 ) is similar to the density used 
bv lMaret et al.l ()2004l ) for reproducing the H2CO emis- 
sion with a non-LTE LVG analysis. The densities used in 
Cases 2 and 3 are slightly lower than the density in the 
hot corino region (where the te mperature is higher tha n 
100 K) of IRAS2A derived by Uorgensen et al.l (J2002I) . 
Higher densities do not reproduce the observed HDO 
emission (the reduced x 2 increases to values much higher 
than 1). Regardless of the density, the derived column 
density of H^ 8 is equal to 6 - 7 x 10 16 cm" 2 . Note 
that at rtH = 6 x 10 5 cm" 3 , the line weakly masers (see 
also [Ncufcld a nd Melnickl I1991D . The column densities 
we obtain are slightly higher , by a factor of two, than 
that derived bv lPersson et al.l (J20121 ). The difference can, 
therefore, come from a combination of the LTE versus 
non-LTE population, gas temperature and line opacity 
(in our model, it is 1.4). The low temperature could in- 
dicate that the gas is thermally decoupled from the dust. 
V(H2 6 0) can then be derived by assuming an isotopic 
abundance ratio | 6 Q/ 18 Q of 560 dWilson et al.lll994D and 
an opr of 3 (see lEmprechtinger et al.l 120101 12013). De- 
pending on «h, we derive an [HDO]/[H 2 0] abundance 
ratio between 0.3 and 8 % (see Table [3]). 

IRAS4A— For IRAS4A, no other HDO lines but the line 
observed in this work are available. The flux of the HDO 
and H^O lines are, therefore, compared with the predic- 
tions obtained by using the same set of physical condi- 
tions as for IRAS2A. We also used another set of param- 
eters presenting a larger source size 9 K of 0.8", consis tent 
with the upper limit given by Pcrs son et al.l (|2012f ) for 
the H^O transition. The increase in 9 S slightly decreases 
the column density of HDO and H^O by approximately 
the same factor (2-3 ), giving similar results to those by 
IPersson et all (pOll . The [HDO]/[H 2 0] ratio, therefore, 
decreases by a factor of two, at maximum. For both sets 
of physical conditions, we predict an [HDO]/[H20] abun- 
dance ratio between 0.5 and 3 % (see Table [3]). 

4. DISCUSSION AND CONCLUSIONS 

The first result of this work is the relatively high wa- 
ter deuteration, ~ 1 %, in IRAS2A and IRAS4A. In 
IRAS2A, t his value i s com patible with the lower limit 
derived by ILiu et al.l f|2011h in the same source (see In- 
troduction). In IRAS4A, this is the first published esti- 
mate. 

Second, as in IRAS16293 and L1157-B1, the water 
deuteration is lower, by about one order of magnitude, 
than the deuteration of formaldehyde and methanol in 
the sa me sources, previously measured by IParise et al.l 
(200j). 

Third, the water deuteration in IRAS2A and IRAS4A 
is very similar to that measured in IRAS16293 by 
iCoutens et alJ (12012T). ~ 3 %, b ut higher than the ratio 
derived bv lPersson et al.l (J2013I ) in the same source. The 
difference between the two results m i ght co me from the 
choice of the method. IPersson et al.1 (|2013T ) derived the 
[HDO]/[H20] ratio from few lines by assuming LTE pop- 
ulation and optically thin emission, whereas the quoted 
column density implies a line opacity ~ 5 and the 203 



TABLE 3 

LVG BEST FIT PARAMETERS FOR HDO AND H^O EMISSIONS. 



Model: 20 K Model: 10 K L1157-B1 IRAS 16293 IRAS2A IRAS4A IRAS4B 



Case 

Density (cm -3 ) 



1 
6 x 10 5 



2 
2 x 10 7 



3 

1 x 10 8 



IRAS2A (T ■■ 



AT(HDO) (cm"" 3 ) 
r(HDO 4 2 ,2-4 2 ,3) 
AT(p-H£ 8 0) (cm" 2 ) 
r(p-HfO) 
AT(H 2 0) (cm- 2 ) 
HDO/H 2 Q 



80 K, 
rg — 



1 x 10 
23 
6 x 10 16 

-0.4 

1.3 x 10 20 

0.08 



h = 0-4 ") 
1 x 10 1S 

2 
7 x 10 16 

0.8 

1.6 x 10 20 

0.006 



6 x 10 iY 

1 

7 x 10 16 

1.5 

1.6 x 10 20 

0.003 



IRAS4A 



AT(HDO) (cm" 2 ) 
r(HDO 4a, 2-4 2 , 3 ) 
A^p-fLfO) (cm" 2 ) 
r(p-H»0) 
AT(H 2 0) (cm" 2 ) 
HDO/H2O 

JV(HDO) (cm" 2 ) 
t(HDO 4 2 ,2-4 2 ,3) 
Af(p-Hi 8 0) (cm" 2 ) 
r(p-H 2 8 0) 
N(H 2 0) (cm~ 2 ) 
HDO/H2O 



1.5 x lO 15 " 

3 
2.5 x 10 16 

-0.4 

5 x 10 19 

0.03 



T = 80 K, S = 0.4 " 



3 x 10 17 " 

0.8 
1.5 x 10 16 

0.05 
3 x 10 19 

0.01 



2 x 10 1Y 

0.4 

1.5 x 10 16 

0.15 

3 x 10 19 
0.007 



T = 

5 x 10 1Y 

0.8 
1.5 x 10 16 

-0.2 

3 x 10 19 

0.016 



80 K, 



0.8 



lx lO 17 "" 

0.3 
6 x 10 15 

-0.02 

1.3 x 10 19 

0.008 



6 x 10 lli 

0.2 
5.5 x 10 15 

0.08 

1.3 x 10 19 

0.005 



GHz line may mase r (see above). On the contrary, 
ICoutens et al.l (|2012l ) uses single-dish observations which 
also encompass the cold envelope even though most of the 
lines have E up > 50K, so that the contamination from 
the outer cold envelope is accounted for. 

The ratio is at least one order of magnitude larger 
than the value measured i n IRA S4B, < 6 x 10~ 4 , by 
|j0rgensen fc van Dishoecld (|2010f ). despite the fact that 
this source lies in the same molecular cloud, NGC1333, 
as IRAS2A an d IRAS4A and it is only ~ 15" away from 
IRAS4A (|Sandellet al.lll991l) . To add to this oddity, the 
deuteration of formaldehyde and methanol in IRAS4B is 
very similar to that measured in IRAS2A and IRAS4. 

Figure [H summarizes the situation, with a plot of 
the measured deuteration of water, formaldehyde and 
methanol in the outflow shock L1157-B1 and in the pro- 
tostars IRAS 16293, IRAS2A, IRAS4A and IRAS4B. In 
the same figure, we also show the theoretic al predictions 
by t he gas-grain model GRAINOBLE (jTaquet et al.1 
2013). Briefly, the model follows the multilayer forma- 
tion of deuterated ices with a pseudo time-dependent 
approach. We report the icy [HDO]/[H2 0] ratio com- 
puted at 3 x 10 5 yr (the typical age of prestellar cores) for 
different constant densities and temperatures, Ay = 10 
mag. The H2 opr, which is difficult to constrain ob- 
servationally, is one of the key parameter s in setting 
the [HDO]/[H 2 0] r atio dTaauet et al.ll20Tl . Following 
the value derived by iDislaire et al.l (|2012| ) towards IRAS 
16293, we used a H2 opr of 10 -3 . The comparison 
between the observations and the theoretical predictions 
shows that the [HDO]/[H 2 0] ratio measured in IRAS2A 
and IRAS4A is reproduced for a large range of physical 
conditions: n H ~ 10 3 - 10 5 cm -3 for T = 10 K and 
n H ~ 10 3 - 10 6 cm- 3 for T = 20 K. On the contrary, 
our mode l cannot reproduce the [HD01/[H 2 0] value re- 
ported bv lJeirgensen fc van Disho cck (2013) for densities 
larger than 10 3 cm- 3 . One possible explanation is that 
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Fig. 2. — Deuterium fractionation of water, formaldehyde, 
methanol of simply (top) and doubly (bottom) deuterated species. 
From left to ( with only one except ion right: predictions of GRAIN- 
OBLE model JTaquet. et al.ll2013l) at 20 K (gr ey) and 10 K (b lack), 
and observations towards L1157-B1 (red, from Codclla et al. 2012), 
IRAS16293 (solid blue, from Loinard et al. 2001; Parisc et al. 2002, 
120041 : IVastel et all I20T7J : [Uoutcns et al. 2012, and dashed blue 
from Persson et al. 2013), IRAS2A and IRA S4A (green, from 
this work, for water, and Parise ct al. 2006, for formaldehyde 
and methanol), and IRAS4B (purple, from Parisc ct al. 2006; 
U0rgensen &: vanD ishocck 2010). 



water ice has formed at a lower H 2 opr (see lTaauet et al.l 
120131 ) or the model is missing some ingredients regarding 
the deuterated ice formation. 

As in our previous work, we note that the larger deuter- 
ation of formaldehyde and methanol testifies to a forma- 
tion of these species on the grain surfaces at a later and 
higher density stage than water, likely the prestellar core 
phase. 

Finally, NGC1333 is a very active star forming re- 
gion undergoing the destruction and alteration from var- 
ious outflows of the first generation stars that might 
have initiated the f o rmation of IRAS2A and IRAS4A 
(iLiseau et al I 119881: iWarin et al.l 119961: iLefloch et all 
119981: IKneefc Sandell |2000|) whereas th e cloud con- 
tainin g IRAS16293 is relatively quiescent (|Mizuno et al.l 
1990). Nevertheless, the similar deuterium fractionation 
derived in IRAS2A, IRAS4A and IRAS 16293 suggests 
that these protostars have followed a similar chemical 
history even though they are located in very different en- 
vironments. However, the [HDO]/[H 2 0] ratio observed 



in IRAS4B by Uorgensen fc van Dishoeckl ()201ClO and in 
IRAS 16293 bv JPersson et all (120131 ) does not fit with 
this conclusion and remains puzzling. 
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